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IN THIS ISSUE... 


In Polar Coordinates and the Cardioid Body of 


M (page 4), Ralph Walde shows that the “body” of 
the Mandelbrot set is the classical cardioid figure, 
and its “head” is a perfect circle. 

This seems remarkable, since those features 
seem at first glance to be subjective, defined only 
by the beholder’s eye — but in fact they are objec- 
tive, and can be characterized by the limit cycles of 
their points. Thus: points in the body of ™ are 
those attracted to a single point, and points in the 
head are those attracted to a two-point cycle. 

Stanley S. Miller likes Steel-Wool Tori (page 2) 
— which are not really fractal: they can get pretty 
kinky, but not infinitely kinky. Miller’s tori are 
quite classical, recalling stately 19th century pen- 
dula dribbling sand on the floor. 


THE SLIDES (S21) 


All four slides accompanying this issue were 
created by Lee Skinner using FRACTINT (version 
13). Three of them have generating functions other 


than the usual f(z) = z +c. Later in this issue, 
Lee shows how to program these functions ( TRIG- 
ONOMETRIC IDENTITIES, page 8). 


Lee uses an 80386 machine with math chip run- 
ning at 16 MHz. The images were photographed 
from the screen of a Super VGA display having 640 
x 480 pixels and 256 colors. 

All the images have a maximum dwell of 1023. 


4516: PLUMAGE. A dayglow cactus. This image is 
part of a #set generated using the exponential 


function f(z) = Ae*+c, A = 0.379 + 0.4791. 
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Center = 2.11434 + 0.40892 i, magnification = 
9.09, aspect ratio = 1.35. 


4517: BIGTHREE. Part of a -set generated using 
the function f(z) = Acosh z+ c, with A = 1+.77. 
An unusual feature of this image is the 
penetration of one family of contours by 
another. 

Center = -1.5708, mag. = 6.94, aspect = 1.35. 


4524: CRESCENT. This is a low-magnification 
M-set image which has been inverted relative to 
a circle of radius r = 0.5 with center 
Do = —0.505, by transforming each point p of 
the complex plane to an inverse point q. 
Expressing p in polar coordinates relative to the 


center pọ we have p= p+ pe” and 
1 ie 


Center = -1.5, mag. = 0.67, aspect = 1.33. 
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4536: STNGRAYP. This unusual 
generated using Ackermann’s 
a 


image was 
generalized 


exponential function z* , where each turn of 
the iteration “crank” adds one more level of 
exponentiation. The iteration function is 


f&2) = č, where we stop iterating when 


a A. 

Center = -1.6826 + -0.07521, mag. = 1.86, aspect 

= 1.34. 

Please note that we’re now calling the slide sets 
“S21” “ete. (CSlide set 219 instead: OF “C20” sete. 
(“Color slide set 20”). 


SUBMITTING ARTICLES 
— Editor 


Here are some guidelines for submitting articles 
for publication in Amygdala. 


1. Type your article as you would like to see it 
appear in the newsletter. Please do not send hand- 
written drafts! 


In order of preference: 


la. On 3.5” diskette in Macintosh format for 
FrameMaker, MacWrite, WriteNow, or MS Word; 
or as a text file. Please also enclose paper copy 
so I can see your intent. 


1b. On 3.5” diskette in IBM format, text file. I 
have no way to deal with 5.25” diskettes. 


lc. Paper copy. 


2. Illustrations should be either grayscale (suitable 
for halftoning) or black/white; not color! 


2a. Normally, illustrations will be printed in full 
column width, so you should make them 3.25” 
wide, if possible — provided that they’re 300 
dpi resolution. If they're grainier, make them 
larger if possible, so that they'll look good when 
reduced to 300 dpi. 


2b. Make sure that you clearly indicate which 
illustrations go where in the text! 


2c. All in all, it’s better not to have captions 
welded into your pictures. Let me put them in 
ad lib. 
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2d. I can handle illustrations on diskette in Mac- 
Paint, MacDraw, CricketDraw, or Adobe Illustra- 
tor formats. 


3. Please send along a short biographical note, 
which I will try to publish in the same issue as 
your article. 


4. Please include your telephone number, in case 
I have to reach you in a hurry with questions. 


FRACTAL GRAPHICS WORKSHOP 


A workshop, Fractal Graphics as a New Art 
Form, was held on August 27-31 at the Glacier 
Institute in the Taos Ski Valley. Situated in the 
beautiful Sangre de Cristo mountains of northern 
New Mexico, the Glacier Institute offers a refresh- 
ing and cool setting for summer conferences. 

Participants came from seven states around the 
country. Rollo Silver made a great presentation, 
with tutorials and lectures which acted as a spring- 
board for cross-fertilization of new ideas and tech- 
niques. His Cibachrome images on display were 
quite inspirational. 

Everyone thoroughly enjoyed the interaction 
and sharing of ideas with people of similar inter- 
ests yet different experiences. Plans are being 
made for another fractal graphics workshop next 
summer during July. 


Dr. Larry Singer, 
Director, Glacier Institute 


STEEL-WOOL TORUS 
— Stanley S. Miller 


Its amazing what a computer can do to perk up 
a harmless branch of 19th century mathematics. 
Lissajous used a narrow stream of sand pouring 
from the base of a compound pendulum to pro- 
duce his curves; I’ve got an IBM clone with an 
8087 co-processor. 


Primitive fragments of Lissajous or Bowditch 
curves were known before the 19th century, of 
course. That is where they picked up such splen- 
did names as “The Witch of Agnesi” or “The Devil’s 
Curve”. What they are, essentially, are parametric 
equations (XY curves) based on trigonometric 
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functions. If you match x = sinit against y = cost 
and run ftthrough a cycle from 0.01 to 6.3, you will 
find yourself looking at a circle. If you then add a 
second component, for example 
x= sint—0.3cos3t and match it against 
y = —0.1sin7¢, you will get a more interesting fig- 
ure, a circle with wiggles. 


I came across this world when I purchased a 
piece of software called CURVES87 for $35 from 
Bridge Software (POB 118; New Town Branch; 
Boston, MA 02258). I began experimenting with 
variations to the parameters and found that with 
some judicious tinkering I could get whole net- 
works that had never been seen in captivity before. 
Our mutual friend Dewdney published one of 
them in the May, 1988 issue of Scientific American, 
and Science Digest published some more in August 
1988. 


The secret of this new approach to the old art is 
to shift the locus of the moving point t off its cycle 
so that it doesn’t get back to the starting point. I do 
this by setting the first element of the first equation 
to sin1.01¢ and the first element of the second 
equation to cos0.99¢. If I then add a second ele- 
ment based on an odd number such as sin3¢ or 
cos7 t, I can spice up the shape. I can get my wig- 
gles off cycle as well by using  sin3.01¢ or 
cos6.99t. I apply these second (or third) elements 
in moderation by slapping a factor such as 0.1 or 
0.3 onto them. Once I throw the equation off bal- 
ance, however, it has to cycle through more 
rounds, usually 315 cycles (m times 100) or 630 
cycles (2m times 100). Fortunately, the most time 
these cycles take to compute is a couple of min- 
utes: we are not in the twentieth century world of 
the Mandelbrot set with its huge drain of computer 
time, because we are not working pixel by pixel. 


1. STEEL-WOOL TORUS (Figure 1, next col- 
umn): 


x = sin0d.99t—0.9sin1.01t+ 0.6sin0.99t 
— 0.4sin2¢t+ 0.1sin50¢t 
y = cos0.99t— 0.9c0s1.01t+ 0.6c0s0.99t 


— 0.4cos2t+0.1cos50t 


Here is a picture of the steel-wool torus, and 
two more samples of the Miller Set to try out: 
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Figure 1 — STEEL-WOOL TORUS 


2. MADNESS, the curve published by Dewdney: 
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Figure 2 — MADNESS 

x = sin0.99t—0.7 cos3.01¢ 

y = cos1.01¢+ 0.1sin15.03¢ 
Upper left x = -1.7679-5908-98 
Upper left y = 1.2072-7858-51 
Lower right x = -CUpper left x) 
Lower right y = -(Upper left y) 
Center: x = y = 0, 
H_stretch = 90.2170-1968-32898 
V_stretch = 82.4167:6878-19051 
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Parameter increment size = 0.01 


3. LOLLY’S POLYP, pictured in Science Digest. 
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Figure 3 — LOLLY’S POLYP 


x = 1.6sin0.99t-—0.9sin1.01t 
— 0.4sin3t+ 0.1sin21¢ 
y= 1.6: "0S 0,991 = 0.908.101 
— 0.4cos3t+ 0.1 cos21t 
Upper left x = -2.8819-7869-72 
Upper left y = 3.0787-2928:39 
-CUpper left x) 
-(Upper left y) 
Center: x = y =0, 
H_stretch = 55.3439-2053-48688 
V_stretch = 32.3185-2846-60037 
Parameter increment size = 0.015 


Lower right x 


Lower right y 


How many more can you find? I’ve got closets 
full of them. 


POLAR COORDINATES AND THE 
CARDIOID BODY OF m 
Ralph E. Walde 


The article by Frank A. Browne in Amygdala #2 
gave a very nice algebraic derivation of equations 
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defining the boundary of the main region of the 
Mandelbrot set. For readers who are not aware of 


why this boundary is called a cardioid', the fol- 
lowing tutorial on polar coordinates may be help- 
ful. Using the polar coordinate definition of a 
cardioid and the relationship of the Mandelbrot set 
to Julia sets, it will be shown that the boundary of 
the “body” of the Mandelbrot set is precisely a car- 
dioid and that the boundary of the “head” is pre- 
cisely a circle. In particular, a trigonometric 
description of the boundary of the body is the set 


of complex numbers {z}, where the real and 
imaginary parts of 2) are given by: 

RZ = 0.25 + 0.5cosO (1 — cos@) 

S zə = 0.5 sin (1 — C056) 


for 0<0<2r, and the boundary of the head is 
given by: 
Rzy = —1+0.25cos8 


for 0 <0 <2n. 


Figure 1 
Boundary of the Body and Head of M 


Polar Coordinates 


Cardioids are (or used to be) defined in calculus 
texts in terms of polar coordinates. Polar coordi- 
nates are an alternate way of describing points and 
curves in the complex plane. Each point in the 


1. Pronounced car dee oyd 


complex plane is specified by a pair (R, ©) of real 
numbers where R is the distance of the point from 
the origin and @ is the angle in radians between the 
positive real axis and the line segment from the 


point to the origin. (R, O specifies the same point 
as the complex number Rcos0+ iRsinð. 


Thus (R,® = (0,0) is the origin 0, (R,98) = 
(1,7/2) is í (R,8 = (1,7) is-1,and (R,9) = 
(./2,n/4) is 1 +17. 


Curves in the complex plane can be described 
in polar coordinates by specifying a relationship 
between R and @. For example, R = r for a posi- 
tive constant r defines a circle of radius r centered 
at the origin, consisting of all complex numbers 
2 = COSO + ir sing. 


The line ax+ by = c is described in polar 
C 
acosð +bsinð 
note that since x = Rcos® and y = Rsinð, we 

have 


coordinates by œR = To see this, 


c= ax+ by = aRcos®ð + bRsinð 
= R(acos® + bsin®) 


Bie a 
acos® + bsin®@ 


Certain curves are far easier to define in polar 
coordinates than in xy-coordinates. Of greatest 
interest to us here is the cardioid. It is a curve 
defined in polar coordinates by 


R = atl —cosd) 


where a is an arbitrary positive constant, which 
controls the size of the cardioid. 


For example, taking a= 0.5, the cardioid 


defined by R= 0.5(1-—cos®@) consists of the 
complex numbers {Zə} , where 
Nz, = 0.5cos0 (1— cos) [1] 
S za = 0.5sin (1 — cos8) [2] 


for 0<$@< 2r. This particular cardioid is shown in 
Figure 2: 


S07 os0. 540.25 


Figure 2 — Cardioid (R = 0.5 (1 -— cos®@) ) 


The curve called the limaçon (lim' uh sun) is 
defined by R = a+ bcos@ for positive a and b. It 
is a generalization of the cardioid, which is a 
limaçon with a = b. Figure 3 shows a limaçon 
with @ = 0.25 and b = 0.75: 


Figure 3 — Limacon (R = 0.25 + 0.75cos0) 


For any positive constant a, R = aĝ defines a 
spiral, called the spiral of Archimedes (Figure 4): 


Figure 4 — Spiral of Archimedes (R = 0.5 0) 
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and R = acos20 defines a four leaved rose (Fig- 
Wire >: 


Figure 5 — Four-leaved Rose (R = acos20) 


See older calculus texts for additional examples of 
curves defined by equations in polar coordinates. 


Regions of the Mandelbrot set 


In looking at Julia sets associated with the com- 


plex function f(z) = rhe C, it has been found use- 
ful to look at regions consisting of those values of 
the complex number c for which f(z) has attrac- 
tive periodic points of a certain period. 


The Body 


Thus, we can define the body of the Mandelbrot 
set to be the region consisting of those complex 
numbers c for which f(z) has an “attractive fixed 
point”, that is, it consists of those c for which there 
is a complex number z, with /(z,) = z, Cixed 


point) and |f'(z,)| <1 
f'\(2) = 2z is the derivative of f(z). 


The boundary of the body is then just those 
fixed points for which the limiting condition 
|Æ (z| = 1 holds, i.e. the set of those complex 


(attractive), where 


numbers c for which there isa z, with z, = z tó 
anid. =f 2) = Rar fe. (2) = 172... which 
means that z, lies on a circle of radius 1/2 cen- 
tered at the origin, so Zz, = 0.5cos® + 0.5 isin, for 


some real 8 with O<O<2n. 
The corresponding value of c is found using 
F : 2 ; 
Z = z + c, which gives c = Z,-—2Z,, or, taking 
parts separately, 


real and imaginary 


Rc = 0.5cos8 — 0.25cos*0 + 0.25sin?0, so apply- 
ing sin0 = 1—cos’0, we get: 
Rc = 0.25 + 0.5cos9 (1 — cos) . [3] 
The imaginary part is given by: 
Sc = 0.5sinð —0.5sinOcos8, or: 
Sc = 0.5sin8 (1 -— cos@) [4] 


We note that these equations [3] and [4] are 
almost the same as [1] and [2], which define the 
cardioid with a = 1/2 — except that 1/4 has 
been added to Nc. Thus we have shown that as 0 
varies between 0 and 2 7, the fixed point z, travels 
around a circle of radius 1/2 centered at the origin, 
and the complex number con the boundary of the 
body of the Mandelbrot set travels around a car- 
dioid with a = 1/2 that has been translated a dis- 
tance 1/4 to the right on the real axis ( Figure 1, 
right). 


The Head 


In a similar manner we can define the head of 
M to be the region consisting of those complex 
numbers c for which f(z) has two points for its 
attractive cycle. For a point c on the boundary of 
the head there are two points 2, and Zz, with: 


Z=2te [5] 

Zaz +é [6] 

and D f'(z)| = 1, that is (since 
EL) = 22): 

|Z,Z,| = 174 [7] 


To determine the values of cfor which the latter 
situation holds, subtract [5] from [6] to get 


43-2) = Zaz, = (2-2) -(%4+2%) 
Or 
(Zi 2) (eee) = 0 [8] 


Since z and z, are distinct, z} — z, cannot be 
zero, so the other factor of [8] must be: 
Z+Z,+1=0 [9] 


Or: 


ey 
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BO ee es, [10] 
Substituting [10] in [7] we get: 


1/4 = |z,z,| = |z,(-1-2,)| = |- z-z 


Applying [6] and [9]: 
1/4= |= %-Z%4+c) = |1+c). 
Thus the boundary of the head is the set of c for 
which |c- (—1)| = 1/4, which defines a circle of 
radius 1/4 centered at —1 (Figure 1, left). 


AUTHOR 
Ralph E. Walde 


I am a mathematics professor at Trinity College 
in Hartford, Connecticut. I am presently teaching a 
combination of mathematics and computer science 
courses. I became interested in fractals and chaos a 
number of years ago and have enjoyed generating 
a few fractal images on a Macintosh Plus computer 
with programs written in Pascal. I am currently 
interested in experimenting with similar programs 
on the Sun 360 color work stations recently pur- 
chased by Trinity. 


TRIVIAL PURSUIT: “AMYGDALA” 


I don’t know who sent in the following items, 
since the signature was cut off the bottom of the 


copy. 


AMYGDALA... 


...1S a medical term describing a part of the human 
brain. The following is quoted from Dorland’s 
Illustrated Medical Dictionary (26th edition, Saun- 
ders, Philadelphia, 1985): 


Corpus amygdaloideum — the amygdaloid body: a 
small, ovoid complex of nuclei partly covered by 
the pyriform cortex, within the tip of the temporal 
lobe, anterior to the inferior horn of the lateral ven- 
tricle of the brain; it is part of the limbic system, 
and is sometimes classified as one of the basal gan- 
glia. The amygdala has olfactory connections, is 
reciprocally connected to the limbic cortex, and 
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projects fibers to the hippocampus, the septum, 
the thalamus, and especially to the hypothalamus. 
Called also amygdala and nucleus amygdalae. 


AMYGDALIN... 


. is a biochemical substance. Its formula is 
C6Hs5CH(CN)-O-C 12H21O0 jo. It is a cyanogenic gly- 
coside characteristically found in seeds and other 
plant parts of members of the Rosaceae family, e.g. 


almonds. The following excerpt is taken from a 


German catalog of chemical substances: ! 


Amygdalin. Amygdaloside; mandelonitrile-ß-gentio- 
bioside; D-mandelonitrile-8-D-glucosido-6-8-D-glucoside. 
Colle;NOn; mol. wt. 457.42. C 52.51%, H 5.95%, N 
3.06%, O 38.47%. Occurs in seeds of Rosaceae; prin- 
cipally in almonds. Obtained by extracting defatted 
almond meal repeatedly with boiling 95% alcohol. 
Structure and synthesis: Haworth, Wylam, J. Chem. 
Soc. 123, 3120 (1923); Kuhn, Ber. 56, 857 (1923); Hud- 
son, J. Am. Chem. Soc. 46, 483 (1924); Zemplén, Kunz, 
Ber. 57, 1357 (1924). Is split by emulsin: Amygdalase 
in emulsin produces glucose and mandelonitrile gluco- 
side which are converted by prunase in emulsin to form 
benzaldehyde, glucose (2 mols.), and HCN. Hydrolysis 
with hot dilute mineral acids gives same end products. 
Cf. Viehoever, Mack “Biochemistry of Amygdalin,” 
Am. J. Pharm. 107, 397-450 (1935). 

_CN 

C,H;-CH 

~ac 90,.0-CeH, 305 


glucose glucose 


Trihydrate, C.»H2;NO;;.383H.O, orthorhombic columns 
from water, m. 200°. m. about 220° when anhydrous. 
The once melted and solidified substance melts again 


at 125-130°. CIN — 42° (c = 1, anhydrous basis). 
One gram of the trihydrate dissolves in 12 ml. water, 
in 900 ml. alcohol, in 11 ml. boiling alcohol. Very 
soluble in boiling water, insoluble in ether. pH of satd. 
water soln. about 7. 


AMYGDALINE and AMYGDALOID... 
... both mean “resembling an almond”. 


AMYGDALOIDS... 


are almond-shaped stones, used for seals 
with deeply engraved icons, in the middle and late 
Minoic times. 
Perhaps this little “amygdaclopedia” can be 
extended, 


1. Here, I’ve substituted the entry from The Merck 
Index of Chemicals and Drugs (7th edition (Merck 
& Co., Inc., Rahway, N.J. U.S.A. 1960). — Ed. 


TRIGONOMETRIC IDENTITIES 


— Lee Skinner 


The following trig identities are invaluable for 
coding fractals that use complex-valued transcen- 
dental functions. 


ety = æ. cosy+74- e+ siny 
= (coshx+ sinhx) - (cosy+ i- siny) 


sin (x+ iy) = sinx. coshy+ i- cosx- sinhy 
cos(x +iy) = cosx:coshy — i-sinx -sinhy 
sinh (x+ iy) = sinhx- cosy+ i- coshx- siny 
cosh (x+ iy) = coshx- cosy+ i; sinhx- siny 
log (x+ iy) = 


Y 
= 5 log (a? +97) +i- (atan œ +2krT), 
ke O de E2 es 
tan(x+1y) = 


sinh2y 
cos2x+ cosh2y 


sin2x 
cos2x+ cosh2y 


tanh (x+ ly) = 
sinh2x 
cosh2x+ cos2y 


sin2 y 
cosh2 x+ cos2y 


Z = g ez 


Comment by R.S.: 


Lee used these identities in programming the 
images shown in the slides (S21). 


Watch out for multivaluedness when iterating 
functions involving log z, such as 


kasr rose yig 


The “fracture lines” or “slip planes” apparent in 
images based on multivalued functions are artifacts 
of the inconsistent choice of values of the function 
on the two sides of the discontinuity. 


I don’t have a solution for this problem. Georg 
Friedrich Bernhard Riemann might, but he dead. 


Anybody have an idea of how to render Riemann 
surfaces graphically nicely? 


RECIPE 


— Leonard Herzmark 


This is from my wife’s grandmother: Mandel- 
brodt!! 


3 eggs 1 tsp. vanilla 


1 cup sugar 1 lb. walnuts or pecans! ! 


1 cup Mazola or 2 tsp. baking powder 


peanut oil 


3.5 cups flour 1 - 2 cups white raisins 

Beat eggs & sugar. Continue beating while add- 
ing oil, vanilla, nuts, raisins, flour, & baking pow- 
der. Oil baking sheet — flour hands well, roll 
loaves as cylinders 2" in diameter, then press into 
elliptical cylinders with flat bottom, minor axis 
about 1". Bake loaves 0.75 hr at 350 ° F. Slice into 
cakes about 0.75" thick, separate, bake 15 min 
longer. Ess, ess, mein kind. Enjoy! 


Mandelbrodt 


1. I can’t see how it can be Mandelbrodt without 
almonds — but who am I to argue with some- 
body’s grandmother? — Ed. 


CIRCULATION 
As of October 15, 1990 Amygdala has 735 sub- 


scribers, 284 of whom have the supplemental color 
slide subscription. 


RENEWAL 


For 53 of you subscribers out there this is the 
last issue of your current subscription. I urge you 
to use the enclosed form to renew your subscrip- 
tion promptly to avoid missing anything. 


AMYGDALA, Box 219, San Cristobal, NM 87564 
505/586-0197 


See ORDER FORM for subscription information 
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